Thirty years of nutrients and biogeochemistry in the Norwegian, Greenland and Iceland Seas, 1990–2019

This dataset contains biogeochemical samples analyzed by the Plankton Chemistry Laboratory at the Institute of Marine Research (IMR), from the Norwegian, Greenland and Iceland Seas. Number of surveys and stations have varied greatly over the last 3 decades. IMR is conducting one annual Ecosystem Survey in April-May each year, with multiple trawl surveys and net tows, but only CTD water collections are reported here. This month-long exercise also has companion vessels from Iceland and the Faroe Islands surveying their own territorial waters. Three transects are the core of the time-series, visited multiple times each year (Svinøy-NorthWest, Gimsøy-NorthWest, Bjørnøya-West). On each station, the CTD cast is sampled for dissolved inorganic nutrients (nitrate, nitrite, phosphate, silicate) and phytoplankton chlorophyll-a and phaeopigments (ChlA, Phaeo) at predetermined depths. At times, short-term projects have collected samples for Winkler dissolved oxygen titrations (DOW) and particulate organic carbon and nitrogen (POC, PN) determinations. This unique data set has seen limited use over the years but is a great contribution towards global ocean research and climate change investigations.

www.nature.com/scientificdata www.nature.com/scientificdata/ determinations and Winkler oxygen titrations (DOW). This set of quality-controlled data, from 30 years of ocean surveys conducted by the IMR, has never been published in its entirety.

Methods
Sample collection and analysis. Seawater samples are collected from Niskin-type water bottles at predetermined depths (Table 1) triggered by a CTD (Conductivity Temperature Depth) unit mounted on a rosette. Two types of CTD have been used on our survey ships over the years (see Gundersen et al. 2 for details) and total numbers of stations sampled has also varied over the years (Fig. 1). Additionally, an annual, full-scale survey covering major parts of the Norwegian Sea and adjacent waters has been done each year since the beginning of the millennium, using one or two IMR ships with participating research vessels from Iceland and the Faroe Islands covering their territorial waters. Sample depths (m)   50  30  20  10  5   75  50  30  20  10  5   100  75  50  30  20  10  5   125  100  75  50  30  20  10  5   150  125  100  75  50  30  20  10  5   200  150  125  100  75  50  30  20  10  5   250  200  150  125  100  75  50  30  20  10  5   300  250  200  150  125  100 75  50  30  20  10  5   400  300  250  200  150  100 75  50  30  20  10  5   500  400  300  200  150  100 75  50  30  20  10  5   600  500  400  300  200  150 100 75  50  30  10  5   800  600  500  400  300  200 150 100 75  50  30  www.nature.com/scientificdata www.nature.com/scientificdata/ Dissolved inorganic nutrients (nitrite, nitrate, phosphate, silicate). After three rinses, each water sample (20 mL) was collected in a polyethylene vial. Up until the turn of the century (1999-2002) most nutrient samples were analysed in real time onboard the research vessels and without a poison or preservative agent. As the research cruise activity expanded, the number of automated analysers could no longer match the number of IMR ships operating simultaneously (see Gundersen et al. 2 for details) and nutrient samples were added 200 µL chloroform to retard biological activity and stored at +4 °C for analysis in the home laboratory, usually within 1-6 weeks. Prior to analysis, the samples were acclimated to room temperature as the evaporating chloroform was evacuated by vacuum. A number of Automated Analyzer (AA) systems have been in use over the years (Table 2), based on methods first described by Bendschneider & Robinson 3 and Grasshof 4 , with a number of minor adjustments suggested by the manufacturers (Alpkem, Skalar). The AA system measures nitrate, nitrite, phosphate and silicate. Briefly, nitrate in seawater is reduced to nitrite coupled to a diazonium ion and, in the presence of aromatic amines, the resulting blue azo-dye is determined spectrophotometrically at 540 nm.   . Selected data will therefore appear in the ASCII format with numeral flags (ASCII flag), or as 8-byte numerals in the netCDF format (netCDF flag). See Jaccard et al. 12 for more detailed descriptions of the terms used below. The data header will list all the flags listed, but only QC flags 0-5 (ASCII) and flags 48-53 or 48b-53b (netCDF) were used to evaluate our data.  www.nature.com/scientificdata www.nature.com/scientificdata/ The nitrate concentration is corrected for ambient nitrite (same analytical method as for nitrate, but without cadmium reduction) measured concurrently. Phosphate reacts with molybdate at low pH and the resulting phosphomolybdate is reduced with ascorbic acid to a blue complex measured spectrophotometrically at 810 nm. With the new Skalar AA purchased in 2017 ( Table 2) the phosphomolybdate peak is now measured at 880 nm. Silicate (silicic acid) is reacting to molybdate at low pH and the resulting silico-molybdate is reduced by ascorbic acid to a blue dye measured spectrophotometrically at 810 nm.    (Table 1), filtered onto a 25 mm GFF membrane and stored frozen at −20 °C until analysis in the land-based laboratory. Historically, pigment samples were transported home by one of the cruise participants, www.nature.com/scientificdata www.nature.com/scientificdata/ as hand-luggage in a cooler with frozen cooler-elements. These days, pigment samples are brought back in specially designed coolers, with an internal temperature recorder, that is rated for −20 °C for a minimum of 3 days. In the laboratory, the samples are thawed in 90% acetone, and stored at +4 °C overnight before analysis on a Turner Design 10 AU fluorometer. Phaeopigments (Phaeo) are measured separately from ChlA, in a second reading of the sample after adding 3 drops of a weak acid (5% HCl). The fluorometer is standardized regularly using a solid standard with known fluorescence, and in accordance with Holm-Hansen & Riemann 5 and the manufacturer 6 . Up until 2008, drift in the light-source was monitored annually, but from 2009 the solid standard has been recorded every time the fluorometer is used. www.nature.com/scientificdata www.nature.com/scientificdata/

Particulate organic carbon and nitrogen (POC, PN) samples. A standard volume (265 mL) is col-
lected from each depth (Table 1) and filtered onto a pre-combusted 25 mm GFF membrane (+450 °C, min. 4 h). Each sample is stored frozen at −20 °C in a pre-combusted glass tube until analysis in the land-based laboratory. Preparations of samples and analysis of elemental C and N is described in detail in Grasshof et al. 7 . Briefly, the dried filter-samples are fumed in acid (conc. HCl) in a desiccator for 4-12 h, before they are dried again and packed in a tin-foil capsule. Analysis of POC and PN is performed on an elemental analyzer (see Table 2 for details) and in accordance with the manufacturer's recommendations.

Dissolved oxygen samples by Winkler titrations (DOW).
Samples for dissolved oxygen are collected in volume calibrated glass BOD bottles (approximate volume 125 mL) and filled from bottom up using a Tygon-tubing. The sample is let to overflow approximately 3 times the volume and great care is taken to avoid small air bubbles at the inside of the sample bottle. Thiosulfate titrations of dissolved oxygen are still done as first described by Winkler 8 but the method has seen some updates and improvements over the years [9][10][11] . Grasshof et al. 7 is describing in detail the current method of sample collection, pre-treatment, and titrations of Winkler samples. Briefly, dissolved oxygen reacts with an alkaline solution (Reagent 1) to form a manganese-hydroxy-complex. Under alkaline conditions, the Mn-complex is reacting with the iodide solution (Reagent 2) and let to precipitate at the bottom of the sample bottle. The sample is added 10 N sulphuric acid to dissolve the iodide precipitate (pH = 1-2.5) and the yellow iodine is titrated by thiosulfate to a clear solution. The titrant is standardized by a known concentration of potassium iodate (KIO3) as described by Grasshof et al. 7 .
Quality control. All data were quality controlled (QC) by analysts using quality flags 0-5 (Table 3) in accordance with Jaccard et al. 12 and the OceanSITES Data Format Reference Manual (http://www.oceansites.org/docs/ oceansites_data_format_reference_manual.pdf). Only QC flags 1, 2, and 5 are made available for these data records. "Good data" (flag = 1) are data that passed QC, "unexpected data" (flag = 2) are data that appear not to conform to expected value (see Technical Validation below for details) but we have no reason to exclude them; "corrected data" (flag = 5) are obvious errors based on notes from the sample sheet, where analytically correct data has been relocated to another depth (e.g. where the entire nutrient profile has been mislabelled and logged in "up-side-down"). In cases where single samples appear mislabelled, exhibit analytical errors, or appear to fall outside expected QC envelope, such values are discarded as "bad data" (flags = 3 and 4 are dedicated to doubtful and erroneous data, respectively).

Data Records
Any use of this 1990-2019 compilation of data should be accompanied by a cititaion of this paper, including a proper use of doi-reference (https://doi.org/10.21335/NMDC-482758181) and citation of the actual data 13 . The data are available from the Norwegian Marine Data Center (NMDC) home page at IMR (http://metadata.nmdc.no/metadata-api/landingpage/49ce6ca612889a957aed47019f4a49f3). All data (Table 4) are hosted www.nature.com/scientificdata www.nature.com/scientificdata/ on an OPeNDAP server at IMR and available both as netCDF4 and csv files, including a pdf-link with reference to a brief description of metadata and methods. The netCDF file (Table 5) follow the SeaDataNet format for profile data (https://www.seadatanet.org/Standards/Data-Transport-Formats). SeaDataNet (https://www. seadatanet.org/) is the standard publication format for biogeochemical data at IMR and the specifics of the netCDF data file format can be found in chapter 4.5.1 of the current version (v.1.22) at https://archimer.ifremer. fr/doc/00454/56547/. Key features of SeaDataNet include rigorous use of the Natural Environment Research Council (NERC, UK) libraries (Table 4), including measurement techniques (P01) and use of data units (P06). The file format is self-contained for both datasets and variables, and we have chosen variable names most familiar to the ocean science community (Table 4). All data in this paper are considered public domain and, as such, some of them (chlorophyll, nitrate, phosphate, silicate, and dissolved oxygen) are also submitted as a minor part of the global data sets in the Copernicus Marine Data Store (https://data.marine.copernicus.eu/products). You will have to sign up with the Copernicus Marine Data Service to get access to the Global Ocean -Delayed Mode Biogeochemical product 14 for the November 2022 release of data.
Distribution and number of stations visited has changed over the three decades presented here (Fig. 1). The first decade (1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999) (Fig. 1). Distribution of stations sampled within a year (Fig. 2) shows a maximum coverage in May during all three decades. Each month is not always visited each year (n < 10 in Fig. 2) but April-June, the period of the Ecosystem Survey, shows most stations visited in each of the three decades. This data set contains www.nature.com/scientificdata www.nature.com/scientificdata/ dissolved inorganic nutrient measurements (nitrite, nitrate, phosphate and silicate) and pigments (ChlA and Phaeo) from most sample stations Fig. 3). On some cruises, dissolved oxygen measured by Winkler titrations (DOW) and particulate carbon and nitrogen (POC, PN) were filter-collected between 1990 and 2019 (Fig. 3).

technical Validation
Quality control of large scale, long-term data sets are crucial to account for potential mislabeling, potential errors in storage and handling of the samples leading to contamination, and potential anomalies during analysis of the samples. The Plankton Chemistry Laboratory is currently using a QC-flagging system to account for the quality of the data that is produced (Table 3). Data with flags 1, 2 and 5 are made available for use in this publication, while flags 3 and 4 are deemed compromised and are not included. Due to the many different people participating in the sampling program on IMR cruises, we must accept that minor mistakes can be made. Some of these can be corrected (e.g. a mislabeled depths with a value that clearly belongs somewhere else) and are given Flag = 5, while others are beyond correction (e.g. samples that shows signs of contamination or appears to have been stored too long) and are given flags 3 or 4. An important QC flag is number 2, where we label data points that are outside expected value, but for no apparent reason. Unfortunately, data from samples collected prior to 2010 in our time-series used a different QC definition of flag number 2, and these data were excluded. Therefore, flags number 2, 3 and 4 were routinely excluded from data sets during the first two decades (1990-2009).
Our collection of unfiltered seawater samples, added small aliquots of chloroform, is not commonplace in ocean research. However, we have not detected any difference between filtered and unfiltered nutrient samples www.nature.com/scientificdata www.nature.com/scientificdata/ (see Gundersen et al. 2 for details). Filter-collection and storage of frozen particulate samples (ChlA, Phaeo, POC, PN), and onboard Winkler oxygen titrations (DOW), do follow internationally recognized guidelines for time-series measurements (e.g. 11 ). The Plankton Chemistry Laboratory at IMR maintain quality control of precision and accuracy by daily assessments of analytical standard curves of internal standards. Our laboratory find it crucial to maintain contacts with others laboratories through regional intercalibration studies such as QUASIMEME (http://www.quasimeme.org/) and in global intercalibrations such as the International Ocean Carbon Coordination Project (IOCCP) and the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) (https://repository.oceanbestpractices.org/handle/11329/883).
The laboratory is using a specific set of QC-criteria (details below) and outliers will be appropriately flagged ( Table 3). All depth-profiles of nutrients are expected to fall within a "QC-envelope" as a function of season ( Fig. 4a-d). We also expect a certain molar relationship between most of the macronutrients (Fig. 5a-c). The filter-collected pigments ChlA and Phaeo (Fig. 6a,b), the particulates POC and PN (Fig. 7a,b) and DOW concentrations (Fig. 8) also have an expected depth profile, but with a high seasonal variability in surface waters. Therefore, a wide range of values are expected in surface waters, depending on time of the year and nutrient availability. As the macronutrients nitrate, phosphate and silicate disappear, a seasonal increase in phytoplankton biomass (ChlA, Phaeo, POC, PN) is expected in surface waters. Therefore, low nutrient concentrations in surface waters are most often associated with elevated concentrations of ChlA, Phaeo, POC and PN. We also expect a semi-constant internal relationship between the two pigments ( Fig. 6c) as we are for the particulate elements (Fig. 7c). Elevated rates of photosynthesis in spring-summer will cause a net increase in DOW concentrations (Fig. 8). Extracellular byproducts from photosynthesis include dissolved organic nitrogen, which is subject to decay into ammonium and further bacterial nitrification into nitrite. Therefore, as macronutrients are approaching a minimum in spring-summer, we usually observe an increase in nitrite concentrations in surface waters (Fig. 4b). Prior to flagging of a parameter, a potential outlier will have to stand out in two or more of the QC-plots indicated above. Finally, operators performing the QC-assessments are mindful to avoid "data grooming" (i.e. applying a too strict QC-envelope), as we may then miss minor changes in the data set that could appear significant on longer term (e.g. decadal) temporal scale.